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NEW CONSTRAINTS ON THE GALACTIC BAR 
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ABSTRACT 

Previous work has related the Galactic Bar to structure in the local stellar velocity distribution. 
Here we show that the Bar also influences the spatial gradients of the velocity vector via the Oort 
constants. By numerical integration of test-particles we simulate measurements of the Oort C value 
in a gravitational potential including the Galactic Bar. We account for the observed trend that C 
is increasingly negative for stars with higher velocity dispersion. By comparing measurements of C 
with our simulations we improve on previous models of the Bar, estimating that the Bar pattern 
speed is J7b/^o = 1-87 ± 0.04, where J7o is the local circular frequency, and the Bar angle lies within 
20° ^ 4>q ^ 45°. We find that the Galactic Bar affects measurements of the Oort constants A and B 
less than ~ 2 km/s/kpc for the hot stars. 
Subject headings: 



1. INTRODUCTION 

The Galaxy is often modeled as an axisymmetric disk. 
With the ever increasing proper motion and radial ve- 
locity data, it is now apparent that nonaxisymmetric 
effects cannot be negle cted (nonzero Oort constant C, 
lOlling fc Deh^l20M h ereafter O&D; nonzero vertex 
deviation. iFamaev et al.l l 2005 | ; asymmetri es in the local 
velocity distribution of stars. lDehnenfl998f) . It is still not 
clear, however, what the exact nature of the perturbing 
agent(s) in the Solar neighborhood (SN) is(are). Possible 
candidates are spiral density waves, a central Bar, and a 
triaxial halo. 

Due to our position in the Galactic disk, the properties 
of the Milky Way Bar are hard to observe directly. Hence 
its parameters, such as orientation and pattern speed, 
have only been inferre d indirectly from obse r vations of 
the inner Galaxy (e.g., iBlitz & Soergel fl99lL iWeinberel 
fl99l . 

However, the Bar has also been found to affect the 
local velocity distribution of stars. HIPPARCOS data 
revealed more clearly a stream of old disk stars with an 
asymmetric drift of about 45 km/s and a radial velocity 
u < 0, with u and v positive toward the Galactic cen- 
ter and in the direction of Galactic rotation, respectively. 
This agglomeration of stars has been dubbed the "Her- 
cu les" stre a m or the "^ -ano maly " . The numerical work 
of lDehnenl (|1999l I2000D andlux| ()200lD has shown that 
this stream can be explained as the effect of the Milky 
Way Bar if the Sun is placed just outside the 2:1 Outer 
Lindblad Resonance (OLR). Using high-resolutio n spec - 
tra of nearby F and G dwarf stars, iBensbv et alJ (2007) 
have investigated the detailed abundance and age struc- 
ture of the Hercules stream. Since this stream is com- 
posed of stars of different ages and metallicities pertinent 
to both thin and thick disks, they concluded that a dy- 
namical effect, such as the influence of the Bar, is a more 
likely explanation than a dispersed cluster. 
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Assuming the Galactic Bar affects the shape of the 
distribution function of the old stellar population in the 
SN, an additional constraint on the Bar can be provided 
by considering the values derived for the Oort constant 
C. In other words, in addition to relating the dynamical 
influence of the Bar to the local velocity field, C provides 
a link to the gradients of the velocities as well. The study 
of O&D not only measured a non-zero C, implying the 
presence of non-circular motion in the SN, but also found 
that C is more negative for older and redder stars with a 
larger velocity dispersion. This is surprising as a hotter 
stellar population is expected to have averaged properties 
more nea rly axisymmetric and hen ce, a reduced value of 
\C\ (e.g.. lMinchev fe Quillenll2007l : hereafter Paper I). 

It is the aim of this letter to show it is indeed possi- 
ble to explain the observationally deduced trend for C 
(O&D) by modeling the Milky Way as an exponential 
disk perturbed by a central Bar. By performing this 
exercise, we provide additional constraints on the Bar's 
pattern speed. Models for the structure in the bulge of 
our galaxy are difficult to constrain because of the large 
numbers of degrees of freedom in Bar models. Thus fu- 
ture studies of structure in the Galactic Center will bene- 
fit from tighter constraints on the parameters describing 
the Bar, such as its pattern speed and angle with respect 
to the Sun. 

2. THE OORT CONSTANTS 

We can linearize the local velocity field (e.g., Paper 
I) about the LSR and write the mean radial velocity Vd 
and longitudinal proper motion ~p l as functions of the 
Galactic longitude / as 

= = K + Asm(2l) + Ccos(2l) (1) 
d 

-p[ = B + A cos(2Z) - C sin(2Z) 

where d is the average heliocentric distance of stars, A 
and B are the usual Oort constants, and C and K are 
given by 
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TABLE 1 
Simulation parameters used 



Parameter 


Symbol 


Value 


Solar neighborhood radius 


ro 


1 


Circular velocity at rrj 


v 


1 


Radial velocity dispersion 


CTu(ro) 


0.05uo or 0.18i;o 


<T U scale length 


rv 


0.9r 


Disk scale length 


r p 


0.37r 


Bar strength 




-0.012 


Bar size 




0.8r cr 



Here r and <fi are the usual polar coordinates and v$ = 
vq+v, where vq is the circular velocity at the Solar radius, 
ro- In this work we primarily consider a flat rotation 
curve (RC) (see, however hence the derivatives of 
U0 in the above equations are identical to the derivatives 
of v. C describes the radial shear of the velocity field 
and K its divergence. For an axisymmetric Galaxy we 
expect vanishing values for both C and K 3 . Whereas C 
could be derived from both radial velocities and proper 
motions, K can only be measured from radial velocities, 
in which case accurate distances are also needed. 

A problem with using proper motions data has been 
described by O&D. The authors present an effect which 
arises from the longitudinal variations of the mean stel- 
lar parallax caused by intrinsic density inhomogeneities. 
Together with the reflex of the solar motion these varia- 
tions create contributions to the longitudinal proper mo- 
tions which are indistinguishable from the Oort constants 
at < 20% of their amplitude. O&D corrected for the 
"mode-mixing" effect described above, using the latitu- 
dinal proper motions. The resulting C is found to vary 
approximately linearly with both color and asymmetric 
drift (and thus mean age) from C w km/s/kpc for 
blue samples to C « —10 km/s/kpc for late-type stars 
(see Figs. 6 and 9 in O&D). Since C is related to the 
radial streaming of stars, we expect non-axisymmetric 
structure to mainly affect the low-dispersion population 
which would result in the opposite behavior for C. In 
Paper I we showed that spiral structure failed to explain 
the observed trend of C. 

Note that the Oort constants are not constant unless 
they are measured in the SN. Due to nonaxisymmctries 
they may vary with the position in the Galaxy. Thus, the 
Oort constants have often been called the Oort functions. 

3. THE SIMULATIONS 

We perform 2D test-particle simulations of an initially 
axisymmetric exponential galactic disk. To reproduce 
the observed kinematics of the Milky Way, we use disk 
parameters consistent with observations (see Table [T]). 
A detailed description of our disk model and simulation 
technique can be found in Paper I. We are interested in 
the variation of C with color {B — V), and asymmetric 
drift v a . We simulate variations with color by assuming 
that the velocity dispersion increases from blue to red 
star samples. For a cold disk we start with an initial 
radial velocity dispersion o~ u — 0.05i>o whereas for a hot 
disk we use o~ u — 0.18urj. The background axisymmetric 
potential due to the disk and halo has the form $o(r) = 

3 Note, however, that C and K would also be zero in the presence 
of nonaxisymmetric structure if the Sun happened to be located on 
a symmetry axis. 
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Fig. 1. — Simulated u — v distribution for a Bar pattern speed 
r2 b = 1.87f2o, Bar angle </>o = 35° and a sample depth d = ro/40. 
The initial velocity dispersion is cr u = 0.18-uo and the Bar strength 
is eb = —0.012. Contour levels are equally spaced. The clump iden- 
tified with the Hercules stream is clearly discernible in the lower 
left portion of the plot. This figure is in agreement with the ob- 
served fo(u, v) (Dclmen 1998; Fux 2001) and simulated distribution 
functions of lDehnenl [120001) and IFuxl 120011) . 

Vq log(r), corresponding to a flat RC. 

3.1. The Bar potential 

We model the nonaxisymmetric potential perturbation 
due to the Galactic Bar as a pure quadrupole 

f ( rk ) 3 , r > r h 
$b = A b (e b )cos[2(0-fi b i)] x { 7 _\a 



Here Ab(eb) is the Bar's gravitational potential ampli- 
tude, identical to the same name parameter used by 
iDehnenl ((2000); the strength is specified by eb = —a 
from the same paper. The Bar length is r b = 0.8r cr 
with r cr the Bar corotation radius. The pattern speed, 
Qh is kept constant. We grow the Bar by linearly vary- 
ing its amplitude, eb, from zero to its maximum value 
in four Bar rotation periods. We present our results by 
changing and keeping ro fixed. The 2:1 outer Lind- 
blad resonance (OLR) with the Bar is achieved when 
r^b/ilo = 1 + k/2 ps 1.7, where K is the epicyclic fre- 
quency. We examine a region of parameter space for a 
range of pattern speeds placing the SN just outside the 
OLR. For a given pattern speed one could obtain the 
ratio r /r O LR t hrough r n /ro L R = f2 b /^OLR ~ fl b /1.7. 

In contrast tolDehnenl (120001) an d similar to lFuxl (|2001f ) 
and iMuhlbauer fc Dehnenl (|2003l ). we integrate forward 
in time. To allow for phase mixing to complete after 
the Bar reaches its maximum strength, we wait for 10 
Bar rotations before we start recording the position and 
velocity vectors. In order to improve statistics, positions 
and velocities are time averaged for 10 Bar periods. After 
utilizing the two-fold symmetry of our galaxy we end 
up with ~ 5 x 10 5 particles in a given simulated solar 
neighborhood with maximum radius d max = ro/10. 

4. RESULTS 

First we show that we can reproduce the results of 
Dehnen and Fux for the Hercules stream. We present a 
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Fig. 2. — Each panel shows the variation of the Oort constant C with Bar angle </>o, for a simulation with the parameters given in Table 
[T]and a particular Bar pattern speed, Q^, and a mean sample depth, d. Solid and dotted lines correspond to cold- and hot-disk values, 
respectively. Columns from left to right, show an increasing Q\> in units of Q(j. Note that the OLR is at f2oLR ~ 1-7. Different rows present 
results from samples with different mean heliocentric distance, d. Good matches to the observed trend in C (vanishing value for the cold 
disk and a large negative for the hot one) are achieved for 20° ^ </>o ^ 45° and 1.83 Q^/Qq ^ 1.91. 



simulated u — v velocity distribution for = 1.87f2oj 
<fio = 35°, and a sample depth d = rn/40 in fig. [TJ The 
initial velocity dispersion is cr u — 0.18i>o- This plot is in- 
deed in very good agreem ent with prev ious test-p article 
jDdmcn 2000; lFull200lf) and N-body (|Fuxll200lh simu- 
lations. 

By a quantitative comp arison of the o bserved with 
the simulated distributions, iDehnenl (|2000l ) deduced the 
Milky Way Bar pattern speed to be f2b/^o = 1.85±0.15. 
In this calculation only the local velocity field, i.e., the 
observed u—v distribution, was taken into account. Now, 
if one could also relate the dynamical effect of the Bar to 
the derivatives of the velocities, an additional constraint 
on Bar parameters would be provided. Since the gradi- 
ents of u and v are hard to measure directly, one needs 
an indirect way of achieving this task. The obvious can- 
didates, as anticipated, are the Oort constants C and 
K. Oort's K, however, is hard to measure as mentioned 
above; hence only C, as estimated by O&D, will be em- 
ployed here. By Fourier expansion of eqs. [1] (see Paper 
I for details) we estimate C from our numerically simu- 
lated velocity distributions for the initially cold and hot 
disks. 

4.1. Variation of C with Bar pattern speed and 
orientation 

In Fig. 13.11 we present our results for C as a function 
of the Bar angle, (j>o (the angle by which the Sun's az- 
imuth lags the Bar's major axis). Each column shows 
a simulation with a different pattern speed, indicated in 
each plot. Rows from top to bottom show C as calcu- 
lated from samples at average heliocentric distances cor- 
responding to d — 200, d = 400, and d = 600 pc, for a 
Solar radius ro = 7.8 kpc. Solid and dotted lines repre- 
sent the results for cold and hot disks, respectively. The 
dashed lines indicates C = 0. C is presented in units of 
= v a / r a ■ To make the discussion less cumbersome, we 
write Ch and C c to refer to the values for C as estimated 
from the hot and cold disks, respectively. 



Ch (dotted lines in Fig. 13. ip varies with galactic az- 
imuth as Ch(0o) ~ sin(20o) for a U of the f2b values con- 
sidered. On the other hand, the cold disk values (solid 
line) exhibit different variations, depending on J7b or 
equivalently, on the ratio r^/roLR- Closer to the OLR 
(left columns of Fig. 13. ip . C c (4>o) approaches the func- 
tional behavior of Ch(0o)- Away from the OLR C c (<po) 
is shifted by 90° compared to Ch((/>o). While both cold 
and hot disks yield in increase in the amplitude of C(0o) 
as the pattern speed nears the OLR, the effect on the 
C c (^>o) is much stronger. This is consistent with our 
expectation that the cold disk is affected more by the 
Bar, especially near the OLR. While close to the OLR 
|Ch(</>o)| < |C c (<^o)|, we observe the opposite behavior 
away from it. This could be explained by the results of 
iMuhlbauer fe Dehnenl (|2003r ). where they find that high 
velocity dispersion stars tend to shift the "effective reso- 
nance" radially outwards. 

According to O&D's "mode-mixing" corrected value, 
C «s for the cold sample and decreases to about -10 
km/s for the hot population. Hence we look for loca- 
tions in Fig. 13.11 complying with this requirement. In 
the transition region, 1.83 ^ fib/rio ^ 1.91, where the 
function C c (0o) transforms into its negative, there are 
specific angles which provide good matches for the ob- 
servations, i.e., C c (0o) ~ while Ch(0o) is significantly 
negative. We performed simulations with different pat- 
tern speeds in the range 1.7 ^ ttb/flo ^ 2.5. It was 
found that only in the transition region in Fig. 13. II could 
one achieve a satisfactory match to the observations. Al- 
though the left- and right-most columns are nearly con- 
sistent with our requirement on C, we reject these values 
of £!b for the following reason. O&D estimated a large 
negative C from both the reddest main-sequence stars, 
which are nearby, and the more distant red giants. This 
is inconsistent with the largest mean sample depth shown 
in Fig. with tt h = 1.81,1.93. We conclude that the 
Bar pattern speed must lie in the range 

fi b /°o = 1-87 ±0.04. (5) 
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The Bar angle with respect to the Sun has been pro- 
posed, as derived from IR photometry, to lie in the range 
15° — 4 5° with t he Su n lagging the Bar major axis, 
whereas iDehnenl (|2000t ) found that the Bar reproduced 



the u-anomaly for 10° ^ <^>o ^ 70°. Examining Fig. 13.1 
we find that our requirements on C result in constraining 
the Bar angle in the range 20° ^ 4>o ^ 45°, depending 
on f2b- This result for the dependence of pattern speed 
on Bar ang l e is in a very good agreement with Fig. 10 
by IDehnenl ({20001 ) which shows a linear increase of the 
deriv ed f2b/^q as a function of cf>o- This is remarkable 
since IDehnenl (2000) obtained his results in a completely 
different way than our method here. 

4.2. The effect of the Bar on A and B 

The Oort constants A and B were found to also be 
affected by the Bar, although to a much smaller ex- 
tent. One could estimate the Bar-induced errors as 
A A = A - A axi and AB = B - _B axi . Here A axi , B axi are 
the resulting axisymmctric values when simulations ex- 
cluding the Bar perturbation are performed. Note that 
we require a different simulation with each of the initial 
radial velocity dispersions used in order to estimate the 
asymmetric drift induced errors. We found that for the 
cold disk the Bar causes AA = 3.5±0.7 and AB = 3±0.6 
km/s/kpc, while the hot disk yielded AA — 0.5±1.5 and 
AB = ± 1 km/s/kpc. These values are calculated for 
ro = 7.8 kpc and vo = 220 km/s. Despite the Bar's effect 
on C, the hot population does provide better measure- 
ments for A and B. 

O&D found an approximately linear variation of A and 
B with color with A increasing and B decreasing. The 
Bar induced errors presented here match the observed 
trend of B but have the opposite one for A. It is possible 
that local spiral structure is responsible for this discrep- 
ancy. 

5. CONCLUSION 

We have shown that the Galactic Bar can account for 
a trend seen in measurements of Oort's C, namely a 
more negative C-value with increasing velocity disper- 
sion. This trend is difficult to explain with other dy- 
namical models (e.g., Paper I). By requiring that the 
Bar model accounts for this trend in C, we have im- 
proved the measurement of the Bar pattern speed, find- 
ing £!b/^o = 1-87 ± 0.04. In addition, the Bar angle is 
found to be in the range 20° ^ 0o ^ 45°. Our result 



for f2b lies well within the es t imate by IDehnenl (|2000f ) 
and that bv iDebattista et all (|2002l ). based on OH/IR 
star kinematics. This study provides an improvement on 
the measurement of the Bar pat tern speed by a factor of 
~ 4 compared to previous work (|Dehnenll2000l ). The im- 
proved constraints on Bar parameters should be tested 
by, and will benefit, future studies of the Bar structure 
in the Galactic center region, that will become possi- 
ble with future radial velocity a nd proper motion studies 
(e.g. , BRAVA, iRich et al.ll200l . 

In addition to a flat RC we have also considered a 
power law initial tangential velocity = vo(r/ro)^ with 
(3 = 0.1, —0.1 corresponding to a rising and a declining 
RC, respectively. We found that for (3 = 0.1 no addi- 
tional error in is introduced. However, for (3 = —0.1 
we estimated J7b/^o = 1.85 ±0.06. From Fig. I3.1l we see 
that the main source of error in the Bar pattern speed is 
du e to the uncertain t y in t he Bar angle. A recent work 
by iRattenburv et all (|2007f ) used OGLE-II microlensing 
observations of red clump giants in the Galactic bulge, 
to estimate a bar angle of 24° — 27°. Using this as an 
additional constraint we obtain Ob/^o = 1-84 ± 0.01. 

While we account for the trend of increasingly nega- 
tive C with increasing velocity dispersion, we fail to re- 
produce the size of the C-value measured by O&D. The 
most negative value we predict is about -6 while they 
measure -10 km/s/kpc. We discuss possible reasons for 
this discrepancy: (i) we expect that the magnitude of 
C is related to the fraction of SN stars composing the 
Hercules stream. It is possible that during its formation 
the Bar changed its pattern speed forcing more stars to 
be trapped in the 2:1 OLR, thus increasing this frac- 
tion. This would give rise to a more negative C from hot 
stars, while leaving the cold population unaffected; (ii) 
the numerical model here uses a distance limited sample. 
Modeling a stellar population with a magnitude limited 
sample is more appropriate to compare to the observed 
measurements, (iii) O&D's "mode-mixing" correction is 
only valid if the Sun's motion in the z-direction were as- 
sumed constant. This is invalidated if, for example, a 
local Galactic warp were present. 
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ABSTRACT 

Previous work has related the Galactic Bar to structure in the local stellar velocity distribution. 
Here we show that the Bar also influences the spatial gradients of the velocity vector via the Oort 
constants. By numerical integration of test-particles we simulate measurements of the Oort C value 
in a gravitational potential including the Galactic Bar. We account for the observed trend that C 
is increasingly negative for stars with higher velocity dispersion. By comparing measurements of C 
with our simulations we improve on previous models of the Bar, estimating that the Bar pattern 
speed is fib/^o = 1.87 ± 0.04, where Q,q is the local circular frequency, and the Bar angle lies within 
20° ^ 0o ^ 45°. We find that the Galactic Bar affects measurements of the Oort constants A and B 
less than ~ 2 km/s/kpc for the hot stars. 
Subject headings: 



1. INTRODUCTION 

The Galaxy is often modeled as an axisymmetric disk. 
With the ever increasing proper motion and radial ve- 
locity data, it is now apparent that nonaxisymmctric 
effects cannot be neglected (nonzero Oort constant C, 
Oiling & Dehncn 2003, hereafter O&D; nonzero vertex 
deviation, Famaey et al. 2005; asymmetries in the local 
velocity distribution of stars, Dehnen 1998). It is still not 
clear, however, what the exact nature of the perturbing 
agent(s) in the Solar neighborhood (SN) is(are). Possible 
candidates are spiral density waves, a central Bar, and a 
triaxial halo. 

Due to our position in the Galactic disk, the properties 
of the Milky Way Bar are hard to observe directly. Hence 
its parameters, such as orientation and pattern speed, 
have only been inferred indirectly from observations of 
the inner Galaxy (e.g., Blitz & Spergel 1991; Weinberg 
1992). 

However, the Bar has also been found to affect the 
local velocity distribution of stars. HIPPARCOS data 
revealed more clearly a stream of old disk stars with an 
asymmetric drift of about 45 km/s and a radial velocity 
u < 0, with u and v positive toward the Galactic cen- 
ter and in the direction of Galactic rotation, respectively. 
This agglomeration of stars has been dubbed the "Her- 
cules" stream or the "u-anomaly" . The numerical work 
of Dehncn (1999, 2000) and Fux (2001) has shown that 
this stream can be explained as the effect of the Milky 
Way Bar if the Sun is placed just outside the 2:1 Outer 
Lindblad Resonance (OLR). Using high- resolution spec- 
tra of nearby F and G dwarf stars, Bensby et al. (2007) 
have investigated the detailed abundance and age struc- 
ture of the Hercules stream. Since this stream is com- 
posed of stars of different ages and metallicities pertinent 
to both thin and thick disks, they concluded that a dy- 
namical effect, such as the influence of the Bar, is a more 
likely explanation than a dispersed cluster. 
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Assuming the Galactic Bar affects the shape of the 
distribution function of the old stellar population in the 
SN, an additional constraint on the Bar can be provided 
by considering the values derived for the Oort constant 
C. In other words, in addition to relating the dynamical 
influence of the Bar to the local velocity field, C provides 
a link to the gradients of the velocities as well. The study 
of O&D not only measured a non-zero C, implying the 
presence of non-circular motion in the SN, but also found 
that C is more negative for older and redder stars with a 
larger velocity dispersion. This is surprising as a hotter 
stellar population is expected to have averaged properties 
more nearly axisymmetric and hence, a reduced value of 
\C\ (e.g., Minchev & Quillen 2007; hereafter Paper I). 

It is the aim of this letter to show it is indeed possi- 
ble to explain the observationally deduced trend for C 
(O&D) by modeling the Milky Way as an exponential 
disk perturbed by a central Bar. By performing this 
exercise, we provide additional constraints on the Bar's 
pattern speed. Models for the structure in the bulge of 
our galaxy are difficult to constrain because of the large 
numbers of degrees of freedom in Bar models. Thus fu- 
ture studies of structure in the Galactic Center will bene- 
fit from tighter constraints on the parameters describing 
the Bar, such as its pattern speed and angle with respect 
to the Sun. 

2. THE OORT CONSTANTS 

We can linearize the local velocity field (e.g., Paper 
I) about the LSR and write the mean radial velocity v,i 
and longitudinal proper motion Jl l as functions of the 
Galactic longitude I as 

= =K + Asin(2Z) + Ccos(2Z) (1) 
d 

~p t =B + A cos(2Z) - C sin(2Z) 

where d is the average heliocentric distance of stars, A 
and B are the usual Oort constants, and C and K are 
given by 



2(7=-- + — - 

r dr r d<j> 

r or r o<p 
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TABLE 1 
Simulation parameters used 



Parameter 


Symbol 


Value 


Solar neighborhood radius 


ro 


1 


Circular velocity at rrj 


vo 


1 


Radial velocity dispersion 


CTu(ro) 


0.05v or 0.18vo 


<Tu scale length 


r CT 


0.9r 


Disk scale length 




0.37r 


Bar strength 


£b 


-0.012 


Bar size 


r\, 


0.8r cr 



Here r and 4> are the usual polar coordinates and = 
Vo+v, where vq is the circular velocity at the Solar radius, 
ro. In this work we primarily consider a flat rotation 
curve (RC) (see, however §5), hence the derivatives of 
in the above equations are identical to the derivatives 
of v. C describes the radial shear of the velocity field 
and K its divergence. For an axisymmetric Galaxy we 
expect vanishing values for both C and K 3 . Whereas C 
could be derived from both radial velocities and proper 
motions, K can only be measured from radial velocities, 
in which case accurate distances are also needed. 

A problem with using proper motions data has been 
described by O&D. The authors present an effect which 
arises from the longitudinal variations of the mean stel- 
lar parallax caused by intrinsic density inhomogeneities. 
Together with the reflex of the solar motion these varia- 
tions create contributions to the longitudinal proper mo- 
tions which are indistinguishable from the Oort constants 
at < 20% of their amplitude. O&D corrected for the 
"mode-mixing" effect described above, using the latitu- 
dinal proper motions. The resulting C is found to vary 
approximately linearly with both color and asymmetric 
drift (and thus mean age) from C « km/s/kpc for 
blue samples to C w — 10 km/s/kpc for late- type stars 
(see Figs. 6 and 9 in O&D). Since C is related to the 
radial streaming of stars, we expect non- axisymmetric 
structure to mainly affect the low-dispersion population 
which would result in the opposite behavior for C. In 
Paper I we showed that spiral structure failed to explain 
the observed trend of C. 

Note that the Oort constants are not constant unless 
they are measured in the SN. Due to nonaxisymmetries 
they may vary with the position in the Galaxy. Thus, the 
Oort constants have often been called the Oort functions. 

3. THE SIMULATIONS 

We perform 2D test-particle simulations of an initially 
axisymmetric exponential galactic disk. To reproduce 
the observed kinematics of the Milky Way, we use disk 
parameters consistent with observations (see Table 1). 
A detailed description of our disk model and simulation 
technique can be found in Paper I. We are interested in 
the variation of C with color {B — V), and asymmetric 
drift v a . We simulate variations with color by assuming 
that the velocity dispersion increases from blue to red 
star samples. For a cold disk we start with an initial 
radial velocity dispersion a u = 0.05wo whereas for a hot 
disk we use a u = 0.18vq- The background axisymmetric 

3 Note, however, that C and K would also be zero in the presence 
of nonaxisymmetric structure if the Sun happened to be located on 
a symmetry axis. 
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Fig. 1. — Simulated u — v distribution for a Bar pattern speed 
Sib = 1.87f2rj, Bar angle </>o = 35° and a sample depth d = rrj/40. 
The initial velocity dispersion is <r u = O.lSvo and the Bar strength 
is eb = —0.012. Contour levels are equally spaced. The clump iden- 
tified with the Hercules stream is clearly discernible in the lower 
left portion of the plot. This figure is in agreement with the ob- 
served fo{u, v) (Dchncn 1998; Fux 2001) and simulated distribution 
functions of Dehnen (2000) and Fux (2001). 

potential due to the disk and halo has the form <j>o(r) — 
Uq log(r), corresponding to a flat RC. 

3.1. The Bar potential 

We model the nonaxisymmetric potential perturbation 
due to the Galactic Bar as a pure quadrupole 

f f^f , r > r h 
$b=A b (e b )cos[2(0-n b i)]x ^ lr As 



Here A b (fb) is the Bar's gravitational potential ampli- 
tude, identical to the same name parameter used by 
Dehnen (2000); the strength is specified by eb = —a 
from the same paper. The Bar length is rb = 0.8r cr 
with r cr the Bar corotation radius. The pattern speed, 
fib is kept constant. We grow the Bar by linearly vary- 
ing its amplitude, eb, from zero to its maximum value 
in four Bar rotation periods. We present our results by 
changing Jib and keeping ro fixed. The 2:1 outer Lind- 
blad resonance (OLR) with the Bar is achieved when 
f2b/^o = 1 + k/2 w 1.7, where k is the epicyclic fre- 
quency. We examine a region of parameter space for a 
range of pattern speeds placing the SN just outside the 
OLR. For a given pattern speed one could obtain the 
ratio ro/roLR through r /r O LR = Slb/^OLR « fib/1-7. 

In contrast to Dehnen (2000) and similar to Fux (2001) 
and Miihlbauer & Dehnen (2003), we integrate forward 
in time. To allow for phase mixing to complete after 
the Bar reaches its maximum strength, we wait for 10 
Bar rotations before we start recording the position and 
velocity vectors. In order to improve statistics, positions 
and velocities are time averaged for 10 Bar periods. After 
utilizing the two-fold symmetry of our galaxy we end 
up with ~ 5 x 10 5 particles in a given simulated solar 
neighborhood with maximum radius d max = ro/10. 

4. RESULTS 

First we show that we can reproduce the results of 
Dehnen and Fux for the Hercules stream. We present a 
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Fig. 2. — Each panel shows the variation of the Oort constant C with Bar angle (j>o, for a simulation with the parameters given in Table 
1 and a particular Bar pattern speed, Qb; ancl a mean sample depth, d. Solid and dotted lines correspond to cold- and hot-disk values, 
respectively. Columns from left to right, show an increasing Q^, in units of Qq. Note that the OLR is at f2oLR ~ 1-7. Different rows present 
results from samples with different mean heliocentric distance, d. Good matches to the observed trend in C (vanishing value for the cold 
disk and a large negative for the hot one) are achieved for 20° ^ cpo ^ 45° and 1.83 ^ O b /S7g ^ 1.91. 



simulated u — v velocity distribution for fib = 1.87fio, 

00 = 35°, and a sample depth d = r /40 in fig. 1. The 
initial velocity dispersion is a u = O.lSt'o. This plot is in- 
deed in very good agreement with previous test-particle 
(Dehnen 2000; Fux 2001) and N-body (Fux 2001) simu- 
lations. 

By a quantitative comparison of the observed with 
the simulated distributions, Dehnen (2000) deduced the 
Milky Way Bar pattern speed to be fib/fio = 1.85±0.15. 
In this calculation only the local velocity field, i.e., the 
observed u—v distribution, was taken into account. Now, 
if one could also relate the dynamical effect of the Bar to 
the derivatives of the velocities, an additional constraint 
on Bar parameters would be provided. Since the gradi- 
ents of u and v are hard to measure directly, one needs 
an indirect way of achieving this task. The obvious can- 
didates, as anticipated, are the Oort constants C and 
K . Oort's K, however, is hard to measure as mentioned 
above; hence only C, as estimated by O&D, will be em- 
ployed here. By Fourier expansion of cqs. 1 (see Paper 

1 for details) we estimate C from our numerically simu- 
lated velocity distributions for the initially cold and hot 
disks. 

4.1. Variation of C with Bar pattern speed and 
orientation 

In Fig. 3.1 we present our results for C as a function 
of the Bar angle, <j>o (the angle by which the Sun's az- 
imuth lags the Bar's major axis). Each column shows 
a simulation with a different pattern speed, indicated in 
each plot. Rows from top to bottom show C as calcu- 
lated from samples at average heliocentric distances cor- 
responding to d = 200, d — 400, and d — 600 pc, for a 
Solar radius ro = 7.8 kpc. Solid and dotted lines repre- 
sent the results for cold and hot disks, respectively. The 
dashed lines indicates C = 0. C is presented in units of 
= u o/ r o- To make the discussion less cumbersome, we 
write Ch and C c to refer to the values for C as estimated 
from the hot and cold disks, respectively. 



Ch (dotted lines in Fig. 3.1) varies with galactic az- 
imuth as Ch(0o) ~ sin(20o) for all of the fib values con- 
sidered. On the other hand, the cold disk values (solid 
line) exhibit different variations, depending on fib or 
equivalently, on the ratio r /roLR- Closer to the OLR 
(left columns of Fig. 3.1), C c (4>o) approaches the func- 
tional behavior of Ch(0o)- Away from the OLR C c (</>o) 
is shifted by 90° compared to Ch((fio)- While both cold 
and hot disks yield in increase in the amplitude of C(0o) 
as the pattern speed nears the OLR, the effect on the 
C c (</>o) is much stronger. This is consistent with our 
expectation that the cold disk is affected more by the 
Bar, especially near the OLR. While close to the OLR 
|Ch(</>o)| < | C c ( 0o ) | , we observe the opposite behavior 
away from it. This could be explained by the results of 
Muhlbauer & Dehnen (2003), where they find that high 
velocity dispersion stars tend to shift the "effective reso- 
nance" radially outwards. 

According to O&D's "mode-mixing" corrected value, 
C w for the cold sample and decreases to about -10 
km/s for the hot population. Hence we look for loca- 
tions in Fig. 3.1 complying with this requirement. In 
the transition region, 1.83 ^ fib/fio ^ 1.91, where the 
function C c (<j>o) transforms into its negative, there are 
specific angles which provide good matches for the ob- 
servations, i.e., C c (0o) ~ while Ch(0o) is significantly 
negative. We performed simulations with different pat- 
tern speeds in the range 1.7 ^ fifc/fio ^ 2.5. It was 
found that only in the transition region in Fig. 3.1 could 
one achieve a satisfactory match to the observations. Al- 
though the left- and right-most columns are nearly con- 
sistent with our requirement on C, we reject these values 
of fib for the following reason. O&D estimated a large 
negative C from both the reddest main-sequence stars, 
which are nearby, and the more distant red giants. This 
is inconsistent with the largest mean sample depth shown 
in Fig. 3.1 with fi b = 1.81, 1.93. We conclude that the 
Bar pattern speed must lie in the range 

fi b /fio = 1.87 ±0.04. (5) 
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The Bar angle with respect to the Sun has been pro- 
posed, as derived from IR photometry, to lie in the range 
15° — 45° with the Sun lagging the Bar major axis, 
whereas Dehnen (2000) found that the Bar reproduced 
the u- anomaly for 10° ^ (f>o 70°. Examining Fig. 3.1, 
we find that our requirements on C result in constraining 
the Bar angle in the range 20° ^ <fio ^ 45°, depending 
on fib- This result for the dependence of pattern speed 
on Bar angle is in a very good agreement with Fig. 10 
by Dehnen (2000) which shows a linear increase of the 
derived fib/fio as a function of <j>Q. This is remarkable 
since Dehnen (2000) obtained his results in a completely 
different way than our method here. 

4.2. The effect of the Bar on A and B 

The Oort constants A and B were found to also be 
affected by the Bar, although to a much smaller ex- 
tent. One could estimate the Bar-induced errors as 
A A = A- A axi and AB = B — _B axi . Here A axi , _B axi are 
the resulting axisymmetric values when simulations ex- 
cluding the Bar perturbation are performed. Note that 
we require a different simulation with each of the initial 
radial velocity dispersions used in order to estimate the 
asymmetric drift induced errors. We found that for the 
cold disk the Bar causes AA = 3.5±0.7 and AB = 3±0.6 
km/s/kpc, while the hot disk yielded AA = 0.5 ±1.5 and 
AB = ± 1 km/s/kpc. These values are calculated for 
r = 7.8 kpc and = 220 km/s. Despite the Bar's effect 
on C, the hot population does provide better measure- 
ments for A and B. 

O&D found an approximately linear variation of A and 
B with color with A increasing and B decreasing. The 
Bar induced errors presented here match the observed 
trend of B but have the opposite one for A. It is possible 
that local spiral structure is responsible for this discrep- 
ancy. 

5. CONCLUSION 

We have shown that the Galactic Bar can account for 
a trend seen in measurements of Oort's C, namely a 
more negative C-value with increasing velocity disper- 
sion. This trend is difficult to explain with other dy- 
namical models (e.g., Paper I). By requiring that the 
Bar model accounts for this trend in C, we have im- 
proved the measurement of the Bar pattern speed, find- 
ing fib/fio = 1-87 ± 0.04. In addition, the Bar angle is 
found to be in the range 20° < 0o ^ 45°. Our result 
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for fib lies well within the estimate by Dehnen (2000) 
and that by Debattista et al. (2002), based on OH/IR 
star kinematics. This study provides an improvement on 
the measurement of the Bar pattern speed by a factor of 
~ 4 compared to previous work (Dehnen 2000). The im- 
proved constraints on Bar parameters should be tested 
by, and will benefit, future studies of the Bar structure 
in the Galactic center region, that will become possi- 
ble with future radial velocity and proper motion studies 
(e.g., BRAVA, Rich et al. 2007). 

In addition to a flat RC we have also considered a 
power law initial tangential velocity = Vo(r/ro)^ with 
f3 = 0.1,-0.1 corresponding to a rising and a declining 
RC, respectively. We found that for /? = 0.1 no addi- 
tional error in fi& is introduced. However, for (5 = —0.1 
we estimated fib/fio = 1.85 ±0.06. From Fig. 3.1 we see 
that the main source of error in the Bar pattern speed is 
due to the uncertainty in the Bar angle. A recent work 
by Rattenbury et al. (2007) used OGLE-II microlensing 
observations of red clump giants in the Galactic bulge, 
to estimate a bar angle of 24° — 27°. Using this as an 
additional constraint we obtain fib/fio = 1-84 ± 0.01. 

While we account for the trend of increasingly nega- 
tive C with increasing velocity dispersion, we fail to re- 
produce the size of the C-value measured by O&D. The 
most negative value we predict is about -6 while they 
measure -10 km/s/kpc. We discuss possible reasons for 
this discrepancy: (i) we expect that the magnitude of 
C is related to the fraction of SN stars composing the 
Hercules stream. It is possible that during its formation 
the Bar changed its pattern speed forcing more stars to 
be trapped in the 2:1 OLR, thus increasing this frac- 
tion. This would give rise to a more negative C from hot 
stars, while leaving the cold population unaffected; (ii) 
the numerical model here uses a distance limited sample. 
Modeling a stellar population with a magnitude limited 
sample is more appropriate to compare to the observed 
measurements, (iii) O&D's "mode-mixing" correction is 
only valid if the Sun's motion in the z-direction were as- 
sumed constant. This is invalidated if, for example, a 
local Galactic warp were present. 
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